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Ile | DNERODUET TON 


Current numerical models have difficulty forecasting 
cases of rapidly developing extratropical maritime cyclones. 
Explosive cyclogenesis has been characterized by Sanders and 
Gyakum (1980) as a central surface pressure fall of at 
least 1 mb/h for 24 h (normalized with respect to 60 degrees 
latitude). These systems present considerable potential 
danger to shipping and fixed ocean platforms. The combined 
effects of high winds, heavy seas and reduced visibility 
developing over short time periods have caused substantial 
property damage, injury and loss of life. The loss of the 
Russian trawler Metrostoy [Mariner's Weather Log, 1979] and 
the sinking of the oil rig Ocean Ranger [Lemoyne, 1982] are 
dramatic examples of the destructive nature of such rapidly 
developing systems. 

Major sources of error in numerical models include data 
gaps or errors, poor initialization, truncation errors and 
inaccurate model physics [Haltiner and Williams, 1980]. 
While all four error sources are potentially important in 
modeling explosive maritime cyclones, the limited amount of 
data available over the ocean certainly is one of the most 
Significant. Satellite systems have helped to fill gaps, 
and the impact of these data on forecasts 1s a major thrust 


of current research. Additionally, the temporal and spatial 


CZ 


scales of these systems are smaller than the scales for 
which some of the models were designed, which results in 
poor resolution of storm features. The comparison of two 
analyses of different resolution can illustrate the impor- 
tance of analyses differences... Finally, the role of convec- 
tive processes in the development of explosive maritime 
cyclones is considerable. Errors in convective cumulus 
parameterization in a model will adversely affect a forecast. 
This study will concentrate on a case of explosive mari- 
time cyclogenesis over the North Pacific during the period 
Smelt o-l6 January 1979. Calland (1983) used analyses by 
the European Center for Medium-range Weather Forecasts 
(ECMWF) based on data collected during the First GARP 
(Global Atmospheric Research Program) Global Experiment 
(FGGE) to i WES ELC AE this system. Other analyses available 
during the FGGE period include the Geophysics Fluid Dynamics 
Laboratory (GFDL) analyses and the Goddard Laboratory for 
Atmospheric Sciences (GLAS) analyses. In this thesis, the 
Slo analysis with its different resolution and initializa- 
tion scheme will be examined. Additionally, two GLAS model 
forecasts will be verified. One set (SAT) of prognoses was 
initialized with all available FGGE data, while the other 
(NOSAT) was initialized with all FGGE data except satellite 
data. The main objectives of this thesis are: 


Oo Diagnostic study of a case of explosive maritime 
cyclogenesis using GLAS analyses; 


Se Diagnostic verification of a global forecast model for 
a case of explosive maritime cyclogenesis; 


se 


o Assessment of the impact of satellite data on the 
model performance. 


The diagnostic approach in this thesis is similar to that 
of Calland (1983). GLAS analyses and model prognoses are 
used to examine mass and circulation (vorticity) budgetsiiaems 
the system. Quasi-Lagrangian diagnostics (QLD) as developed 
by Johnson and Downey (1975) are used to conduct the budget 
Studies. The budget volume is defined by a cylindrical 
coordinate system which is centered on the cyclone surface 
center at each time. The advantage of this technique is 
that it focuses on those processes related directly to 
cyclone development. Verification of the GLAS prognoses is 
achieved by comparing both the basic fields and budget sta- 
tistics to those from the GLAS and ECMWF analyses. 

Appendix A is a discussion of the FGGE level II-b and 
level III-b data. Chapter II describes the results of the 
Synoptic diagnostic studies performed on the GLAS analyses 
and model predictions and compares them with the results of 
Calland (1983). Chapters III and IV explore the resultswes 
the cyclone mass and vorticity budget studies with comparison 
to Calland (1983). Chapter V summarizes the thesis conclusions 


and includes recommendations fot tucure seuadies. 
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Peco nc DiSCUSSION 


A. GENERAL 

Two FGGE analyses and two prognoses are compared for the 
13-16 January 1979 North Pacific storm. The two analyses 
are made by the ECMWF and by the GLAS. Both analyses are 
based on FGGE level II-b data. The two GLAS prognoses are 
from a GLAS analysis which included all available FGGE data 
(SAT) and from a separate analysis that excluded satellite 
data (NOSAT). Intercomparison of the GLAS and ECMWF analyses 
will illustrate the impact of differences due to horizontal 
resolution and different objective analysis schemes in por- 
traying explosive maritime systems. The comparison of the 
‘two GLAS prognoses with each analysis produces a qualitative 
evaluation of model performance with different data sets. 
Finally, intercomparison of the SAT and NOSAT prognoses 
assesses the importance of satellite data in producing a 
better forecast. The evaluation will focus on the surface, 


500 mb and 250 mb features. 


B. ANALYSIS INTERCOMPARISONS 

While both the ECMWF and GLAS analyses are based on FGGE 
level II-b data, substantial differences exist between the 
schemes used to produce the level III-b analyses. The ECMWF 
data assimilation system is a three-dimensional, multivariate, 


optimal interpolation scheme using a nonlinear normal-mode 


ae: 


initialization. A 15-level model with horizontal resolution 
of 1.875 degrees of latitude and longitude is used in the 
update cycle. 

The GLAS objective analysis uses a modified Cressman 
(successive correction) scheme based on a method by Berg- 
thorsson and Doos.(1955). The assimilation/forecast model 
consists of nine vertical levels in sigma coordinates with a 
uniform non-staggered grid (4 degrees latitude by 5 degrees 
longitude). Details of the GLAS FGGE level III-b analysis 
can be found in Appendix A. 

The storm tracks for the GLAS analysis, SAT and NOSAT 
prognoses are depicted in Fig. l. Figures appear in Appen- 
dix C. The reader should consult Calland (1983) to observe 
Similar figures for the ECMWF analysis, as they will not be 
repeated in this thesis. Figure numbers ascribed to the ECMWF 
analysis in the foregoing discussion refer to the figures in 
Calland (1983). A comparison of the central surface pressures 
of the two analyses reveals agreement to within 1 mb with the 
exception of 12 GMT 14 January when the ECMWF analysis is 
4 mb deeper than the GLAS analysis. Storm positions of the 
two analyses are consistent throughout the 72-h period, so 
only the storm track from the GLAS analyses is shown in 
MAES ree lly 

At 00 GMT 12 January atmospheric conditions southeast of 
Japan are favorable for frontogenesis. The GLAS analysis 


sea-level pressure (SLP) (Fig. 2c) contains a broad trough 


16 


eerenaing NOrtnheast-Ssouthwese (vicinity of 35 N, 150 E) 
within a checkerboard pattern of strong cyclone and anticy- 
clone features. The major circulation centers include the 
ppeserian high near 55 N, 100 E (off the chart), a filling 
Teme 9oe mo; at o3 N, 170 EF, a large maritime high (1030 mb) 
centered over the north-central Pacific (vicinity of 30 N, 
178 W), and Typhoon Aliceat 14 N, 138 E. The ECMWF surface 
analysis (Fig. 4c, Calland) indicates similar features. 
aes circulation funnels cold, polar continental air south- 
ward over the northwestern Pacific toward the northward 
flowing warm, moist tropical air, which enhances the polar 
front southeast of Japan. 

By 00 GMT 13 January, a closed surface circulation is 
observed in both analyses (Fig. 3c and Fig. 7c, Calland). 
The positions and intensities of the central surface pres- 
sure in the two analyses are in close agreement. Likewise, 
the thermal fields over the cyclone are similar. Both analyses 
indicate a shallow 500 mb trough west of the surface low 
Peeoeeue 135 —. The ECMWF vorticity field (Fig. 7b, Calland) 
shows a maximum associated with the trough as expected. 
However, the GLAS vorticity field (Fig. 3b) has a smaller 
maximum. Throughout the life of the cyclone, the horizontal 
resolution of the GLAS 500 mb vorticity analyses appear to 
produce excessive smoothing with more symmetrical maxima 
than those seen in the ECMWF analysis. The consequences of 


this difference are explored in the discussion of the vorticity 


IL 


budget results. The two 250 analyses are similar, and both 
indicate a jet streak poleward and somewhat west of the 
surface low. 

At 00 GMT 14 January, there are virtually no differences 
between the analyses in the position and intensity of the 
central surface pressure. The GLAS analysis thickness ridge 
(Fig. 4c) remains stronger than on the ECMWF analysis (Fig. 
10c, Calland). The 500 mb analyses are similar, although 
the GLAS analysis (Fig. 4b) has not resolved the vorticity 
maximum migrating toward the system. The positive vorticity 
advection is well represented on the ECMWF analysis (Fig. 
10b, Calland). At 250 mb, the jet maxima is now positioned 
to the southwest of the surface low in both analyses, which 
creates favorable conditions for strong upper-level diver- 
gence and explosive development. 

By 12 GMT 14 January, a larger difference in the two 
Surface pressure analyses is present. The ECMWF analysis 
(Fig. l6c, Calland) places the central pressure at 965 mb, 
and the GLAS analysis (Fig. 5c) has 969 mb. This is the 
first period of rapid deepening of the system, and the GLAS 
analysis presumably suffers from limited data near the storm 
center and an inability to describe the central feature with 
the coarse resolution. Examination of the level II-b data 
available to the GLAS analysis reveals that there are only 
two surface pressure reports (ship) within 5 degrees lati- 


tude of the storm center. A 500 mb trough just to the west 


das 


of the surface center is similarly represented in the two 
analyses (Figs. 5b and 16b, Calland), which indicates the 
self-development of the storm during the prior 12 h. Cyclonic 
flow in the middle troposphere GLAS vorticity field (Fig. 

Sb ) Indiea tes a Maximum over the surface low, as does the 
ECMWF analysis (Fig. 16b, Calland). Due to coarser resolution, 
the GLAS vorticity maxima are not as large, and are generally 
more circular than those in the ECMWF analysis. The jet 

axis is now well south of the surface low, and the upper-air 
flow is similar in both analyses. 

Between 00 GMT 15 January and 16 January, the cyclone is 
in the occlusion phase of development. The fields of the 
two analyses are similar with the exception of the vorticity 
maxima already described. 

The overall comparison of the GLAS and ECMWF analyses 
shows close agreement. The translation of the storm is very 
Similar throughout the three-day period and, with minor 
exceptions, the intensities of the surface low centers are 
bmnane tire decreased resolution of the GLAS analysis is 
most evident in the vorticity field, where ae maxima are 
either not resolved or oddly shaped. The upper troposphere 
winds are very similar throughout the period in the two 
analyses. In summary, the GLAS analysis performs well on 
this system considering its coarser resolution and less 
complex objective analysis scheme. 

Another comparison of the ECMWF and GLAS analyses for 


January 1979 by Baker (1983) revealed generally close agreement. 


ig 


The positions of the low centers were very similar. The 
major difference in that comparison was that the ECMWF 
analysis depicted a central surface pressure 5 mb lower than 
the GLAS analysis on 19 January. This difference was 


attributed tosthe coarser grid used by the GLAS modem 


C. GLAS PROGNOSES VERIFICATION 

The GLAS prognoses with all available data (SAT) and the 
GLAS prognosis with all except satellite data (NOSAT) are 
verified using the GLAS analyses discussed above. Central 
surface pressure in 12-h increments and the storm tracks 
appear in Fig. 1. It was expected that the SAT prognoses 
would outperform the NOSAT prognoses. As can be seen in 
Fig. 1, this is not necessarily true. The storm track in 
the SAT prognosis during the period of 13-16 January is 
considerably -farther from the analysis track than the NOSAT 
prognosis. The SAT prognosis did forecast storm intensity 
more acorn ee than the NOSAT prognosis in the early 
periods, but after the 00 GMT 15 January the SAT prognoses 
are worse than NOSAT both in position and intensity. 

The SAT prognosis starts with the GLAS analysis at 00 GMT 
13 January (Fig. 3). Noteworthy differences in the initial 
analysis for the NOSAT model are seen in the 250 mb flow 
(Fig. 6a), as the removal of the satellite sounding and 
cloud-drift wind data has altered the intensity of the jet 
stream. The 250 mb jet streak at 140 E on the NOSAT initial 


analysis extends farther eastward than the SAT jet streak 
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fice sa), ald winds Gast Of the incipient region are 10-20 
m/s larger in the NOSAT analysis. 

After only 12 hours, noticeable differences develop among 
the SAT and NOSAT forecasts and the analysis. The surface 
pressure of the cyclone depicted in the SAT prognosis has 
been overdeepened by 3 mb. The storm has tracked more 
northeastward than eastward as on the analysis, and is 130 
Rm away from the analysis center (Fig. 1). The NOSAT prog- 
nosis has likewise overdeepened the system by 4 mb. How- 
ever, the position error of the NOSAT prognosis is less than 
80 km. Consistent with the overdevelopment, the 1000-500 
mb thickness ridge of both prognoses is greater than the 
analysis (Figs Saeresecmand 9C)- 

Pe eget tCrence We Observed ak this time period in 
the 250 mb jet: stream. The NOSAT prognosis 250 mb wind 
field (Fig. 9a) shows the 75 m/s jet streak extending farther 
eastward than on the analysis (Fig. 7a). Consequently, the 
Surface low is already positioned under the front left 
quadrant of the jet streak. The SAT prognosis 250 mb winds 
(Fig. 8a) are 10-20% weaker than the analysis. 

At 00 GMT 14 January, the storm in the SAT prognosis 
(Fig. 10c) has continued to overdevelop and is ncw 6 mb 
deeper than the analysis (Fig. 1). Movement during the 
period is similar to the previous period and the storm has 
a 370 km position error. The NOSAT prognosis (Fig. llc) 
also overdevelops the system with a 14 mb difference in 


central surface pressure, but has only a 120 km position 


Ja Al 


error. The thermal ridges in each prognosis continue to be 
Stronger than the analyses, which is consistent with the more 
advanced stage of development. Consistent with the thermal 
ridges, 500 mb flow above the surface low appears to be more 
cyclonic in the prognoses than in the analyses. However, 

5900 mb absolute Gor mietes differences at this time are 
minimal (Figs. 4b, 10b and 11b). Smaller magnitudes are 
still observed in the SAT jet stream velocity (Fig. 10a); 
however, the difference does not appear to be as great as 

in the previous period. 

In the 12-h period ending 12 GMT 14 January, the three 
storm tracks are generally parallel with similar displace- 
ments (Fig. 1). Overdevelopment continues at a consistent 
rate and the SAT storm is 10 mb deeper than the analysis. 
Development in the NOSAT prognosis slows as the storm enters 
the occlusion phase. Nevertheless, the NOSAT surface cen- 
tral pressure is 16 mb deeper than the analysis. The SAT 
and NOSAT thickness fields and 500 mb flow (Figs. 12 and 13) 
continue to indicate a more advanced stage of development 
than in the analysis stemm eee 

The 500 mb analysis contains a short wave trough trailing 
behind the surface low and a thermal trough behind the short 
wave, which suggests continued Be Oriente of the system 
(Fig. 5b). Both the SAT (Fig. 12b) and NOSAT prognoses 
(Fig. 1l3b) indicate a S00 "mow trough cHrlaree over 


the surface low which signifies self-development 
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during the prior 12 h. Also, the SAT and NOSAT prognoses 
have similar thickness ridges (Figs. 12c and 13c). Each 
prognosis indicates a warm air sector wrapping around the 
system center, which is consistent with the onset of 
occlusion. 

The 250 mb jet stream of the SAT prognosis (Fig. 12a) 
is now much more similar to the analysis (Fig. 5a) than in 
previous periods, particularly in the jet streak southwest 
@eecne storm. The NOSAT prognosis 250 mb flow (Fig. 13a) 
differs somewhat from the analysis (Fig. 5a). Flow is of 
Similar mangitude but is less cyclonic over the surface center 
in the prognosis. As the storm continues to move northeast- 
ward, the influence of the jet on development is diminished 
considerably after this period. 

The period ending at 00 GMT 15 January is noteworthy 
because the observed storm (Fig. 14) continues to develop 
rapidly, while in the SAT and NOSAT prognoses (Figs. 15 and 
16) the storm proceeds into the occlusion phase. Observed 
feeropient iS So intense that the large differences between 
the meenoses and the analyzed surface pressure is diminished 
considerably. The analyzed storm center moves twice as far 
aS in either prognosis and increases the distances between 
HOw centers to 692 km (SAT) and 563 km (NOSAT) (Fig. 1). 

The 500 mb analysis (Fig. 14b) shows greater cyclonic 
Curvature than in the prognoses (Figs. 15b and 16b). A major 


reason for the difference in storm translation between the 
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analyses and the prognoses is due to the position of the 
surface low relative to the jet stream. The analyzed low 
remains centered on the jet axis, while the low positions 
depicted in the prognoses are north of the jet axis, out of 
the baroclinic zone and into the colder air. 

By 12 GMT 15 January all three representations of the 
storm are in the occlusion phase of development. The SAT 
(Fig. 18c) and NOSAT (Fig. 19¢c) storms aneg3s andy? tar 
respectively, deeper than in the analysis (Fig. l7c). The 
central surface pressure in the analysis remains nearly the 
Same as in the preceding period, while the surface pressure 
has dropped ll mb and 5 mb in the SAT (Fig. 18c) and N@Sam 
(Fig. 19c) prognoses. This is unexpected inasmuch as the 
storm in both prognoses had occluded earlier than in the 
analysis. Closed circulations are observed at 500 mb in all 
three representations. The relative stages of occlusion are 
evident in the areal extent of the closed circulation, as 
the center in the analysis (Fig. 17b) is smaller than in 
either of the prognosis 500 mb centers (Figs. 18b anda 
Highly cyclonic flow is even seen to a great extent at 250 
mb in each prognosis (Figs. 18a and 19a). 

By 00 GMT 16 January, the storm in both of the prognoses 
and in the analysis has begun filling. The analyzed storm 
is now over the Alaskan coast, and the biggest difference 
between the analysis and the prognoses is in the thickness 


fields (Figs. 20c, 21@ and 22¢).9Warm aitugaecs not wrae 
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around the analysis storm as in the two prognoses. Apart 
from the position of the surface centers (Fig. 1), no other 
Significant differences are noted for this final analysis 


period. 


D. SUMMARY OF SAT/NOSAT DIFFERENCES 

In the Seely) stages of development, the principal differ- 
ence observed between the two prognoses is the weaker jet 
stream winds in the SAT prognosis. One effect is to trans- 
late the system less accurately (Figs. 8a and 9a). The 
level II-b cloud-drift winds computed by the Meteorological 
Satellite Center (MSC) in Japan were found to contain sub- 
stantial errors due to invalid altitude assignments [McPherson, 
1984]. Upper-level winds were corrected at GLAS prior to 
assimilation of the level III-b data (see Appendix A). 
Nevertheless, errors in the corrected cloud-drift winds were 
large enough to decrease the jet intensity and to cause the 
SAT prognosis to translate the storm more slowly than the 
NOSAT prognosis. A difference plot (NOSAT-SAT) of jet 
stream winds at 00 GMT 13 January is an ligt jaalKe ae ao 
The largest difference between the two prognoses (20 m/s) 
coincides with several cloud-drift wind reports that are 
approximately 25 m/s (50 kt) slower than Sumreunettie cai r— 
craft reports. Once the system moves north of the jet axis 
(about 12 GMT 14 January), the SAT prognosis translates the 


storm at much the same rate as the NOSAT (Fig. 1). This is 
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consistent with the observation that the large differences 
in the 250 mb flow are limited primarily to the jet stream. 

For the majority of the 72-h period, the two prognoses 
develop the storm similarly in terms of intensity. Over- 
development of the storm is forecast during the first 36 5h, 
especially in the NOSAT prognosis. In the latter stages of 
development and occlusion, the positions in the two prognoses 
become closer, but the intensities vary considerably as the 
SAT prognosis deepens appreciably while NOSAT deepening slows 
between 00-12 GMT 15 January. Meanwhile, the analyzed storm 
continues to deepen and translate well to the northeast of 
both progqnosess (tidg-). 

In summary, both prognoses initially overdevelop an 
already rapidly developing system. Because of this initial 
overdevelopment, both prognoses occlude too soon and stall 
over the Aleutian Island chain. The inital impression is 
that the presence or omission of upper tropospheric satellite 
winds does have a considerable effect on model performance. 
The fact that neither prognosis accurately forecasts the storm, 
however, indicates that the presence or absence of satellite 
data is not the sole contributor to the forecast errors. 

The performance of operational forecast models in this 
case was guite different than Phat of the GLAS FGGE prognoses. 
Prognoses from the Fleet Numerical Oceanography Center 
(FNOC) underdeveloped the storm, rather than overdeveloping 


it like the GLAS model. Position errors for the FNOC 
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prognoses in the early periods were as large as 900 km. The 
GLAS FGGE prognoses performed well in penne: toon with the 
mieoc prognoses, particularly in the first 48 h. The impact 
of the comprehensive FGGE data was positive in the fore- 


@eseing of this storm. 
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IIi. MASS BUDGET ANALYSIS 


A. GENERAL 

One of the key observations in the synoptic discussion 
(Chapter II) is the largely zonal mid-tropospheric flow over 
the incipient region in the early stages of development. 
Without a mid-level trough to explain cyclogenesis, other 
forcing mechanisms must be explored. Because of the large 
baroclinity over the incipient region, destabilization of 
the surface layer may explain the early stages of storm 
development. Additionally, evaluation of the horizontal and 
vertical mass circulation can substantiate the qualitative 
observations made in Chapter II, since storm intensity and 
pressure tendencies are directly related to the mass flux 
into and out of the budget volume. 

This chapter quantitatively evaluates the destabilization 
and horizontal and vertical mass transport within the budget 
volume. The generalized quasi-Lagrangian budget equate 
and the QLD mass budget equation are included as Tables [I 
and II. Initially, budget results were analyzed for two 
budget volumes: one with a 6 degree radius (666 km) and one 
with a 10 degree radius (1111 km). The purpose of studying 
the two volumes is to investigate the large-scale storm 
Circulation and mass transport (radius 10) and the more 


intense processes closer to the storm center (radius 6) 
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[Calland, 1983]. The grid-size difference between the GLAS 
and ECMWF analyses is depicted in Fig. 24. The radius 6 
budget results for the GLAS 4 x5 degree grid contain many 
inconsistencies due to the grid size. Consequently, this 
Seay will focus on the radius 10 budget results for both 
the mass and vorticity study. As in the previous chapter, 
four data sets are examined. The GLAS analysis is first 
compared with the ECMWF analysis as described by Calland 
(1983), and then the GLAS analysis is used as a basis for 


comparison of the SAT and NOSAT storms. 


B. POTENTIAL TEMPERATURE AND STABILITY ANALYSIS 

As discussed in Chapter II, the storm development 
occurred within a large-scale flow of polar continental air 
around an intense Siberian high that extended over the warmer 
west North Pacific. These conditions produce a destabiliza- 
tion of the boundary layer over the incipient storm region, 
Mmeereis a key element in Maritime cyclogenesis as described 
by Roman (1981) and Sandgathe (1981). The amount of destabi- 
lization in the surface layer depends on both latent and 
sensible heat flux from the ocean, which in turn is largely 
dependent on air-sea temperature difference. 

Seipeart son Or tae radmus 10 ECMWF analysis area-averaged 
Mecential temperature {Fig. 25b) with that of the GLAS analy- 
Sis (Fig. 25a) reveals that mid-tropospheric temperatures of 
the ECMWF analysis are generally several degrees higher than 


the GLAS analysis. Over the 72-h period, 500 mb temperatures 


29 


of the ECMWF analysis decrease from 3l1l K to 294 K and in 
the GLAS analysis from 309 K to 294 K. The ECMWF analysis 
at 1412 indicates more rapid cooling in the 775 and 600 mb 
layers than in the GLAS analysis, which enhances destabili- 
zation. In the final periods, weaker advection of cold air 
into the surface layer of the GLAS analysis volume contributes 
to its destabilization. While the ECMWF analysis has a 
warmer surface layer in the first 48 h, during occ lusdion maaan 
air is advected into the volume at a greater rate, and by 
1600 the surface layer of the ECMWF analysis is cooler than 
in the GLAS analysis. 

Static stability values similar to those defined by 


Sandgathe (1981) as: 
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will be used within the 1000-500 mb layer. The specific 
volume of the layer is assumed constant with time. The 
potential temperature fields indicate rather uniform structure 
for the lower troposphere through time, so the stability 
index can be simplified as - = 

The static stability trends for the GLAS and ECMWF 
analyses are shown in Fig. 26. Both analyses show a destabi-— 
lization of the layer for 36 h followed by an increase in 


stability after 1412. Because temperature differences between 


the two analyses are consistent throughout the lower troposphere 
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itieeme cariy periods, stability indices vary only slightly. 
The effect of the greater cooling between 500-800 mb in the 
ECMWF analysis (Fig. 25b) at 1412 is reflected in the slightly 
lower stability index. Similarly, the higher stability of 

the ECMWF analysis in the final period is a direct reflection 
of the cold air advection into the surface layer observed in 
Bem. 25D. 

Comparison of the SAT and NOSAT potential temperature 
fields (Fig. 27) with the GLAS analysis indicate several 
important differences. The storm representation in the SAT 
prognosis cools in the mid-troposphere between 1300 and 1312, 
although surface layer temperatures are similar to the GLAS 
Mieysis. The strong cooling at 600 mb (500-700 mb layer) 
is responsible for the rapid decrease in stability for the 
SAT storm at 1312 (Fig. 26). Mid-tropospheric temperatures 
in the SAT and NOSAT prognoses remain lower than in the 
GLAS analysis throughout the storm and the stability like- 
itisc remains low. Low stability values computed for the 
NOSAT prognosis result from a combination of slightly lower 
mid-tropospheric temperatures and slightly higher surface 
layer temperatures as compared with the GLAS analysis. The 
low stability found in the early periods of both prognoses 
May contribute to convective processes in the early periods, 
which would help explain the overdevelopment in the first 
24 h of both prognoses. 

fovestaoumwitvy women the Only facter Contributing to the 


rapid deepening in the prognoses. The NOSAT storm at 00 GMT 


Sl 


14 January is 8 mb deeper than the SAT storm in spite of the 
lower stability for the SAT storm. Throughout the storm, 
both prognoses indicate lower stability (1 K/100 mb) than 

the analyses. In the early periods, this is due primarily 

to lower 500 mb temperatures, which is consistent with the 
respective storm tracks toward the cold air. In the later 
periods, this effect in combination with a warm surface layer 


tends to further destabilize the lower troposphere. 


C. VERTICAL MASS PRANSE CR. 

The storm vertical mass transport is examined through 
vertical time sections of omega fields. For this study, 
omega is calculated using the kinematic method. The O'Brien 
(1970) correction scheme is utilized to correct for residual 
imbalances. 

Time sections of area-averaged omega in the GLAS (Fig. 
28a) and ECMWF analyses (Fig. 28b) each indicate mid- 
tropospheric maxima during the explosive phase of development, 
although there are some differences in magnitude. The 
GLAS analysis indicates a slightly larger vertical velocity 
maximum that lags by 12 h the ECMWF maximum. These time 
differences in maximum vertical motion are consistent with 
the central surface pressure of the ECMWF being 4 mb deeper 
than the GLAS analysis at 12 GMT 14 January. During the 
subsequent 12-h period, when vertical motion of the GLAS 


analysis was at a maximum, a 20 mb pressure fall was analyzed. 


ee 


The omega time-sections for the SAT and NOSAT prognoses 
eeegeeos)  foth indicate larger upward vertical motion during 
the early periods than in the analyses. The NOSAT prognosis 
shows stronger omega over a longer period of time, which is 


consistent with more rapid development (Fig. 1). 


D. HORIZONTAL MASS FLUX 

Vertical time sections of horizontal mass transport for 
the four data sets are depicted in Figs. 30 and 31. fMThe 
GLAS analysis (Fig. 30a) has a two-layer regime of low-level 
mass influx associated with surface layer convergence and a 
layer of mass outflow associated with upper-level divergence. 
The surface layer convergence, while more intense than the 
upper-level divergence, is confined to a more shallow layer. 
The weaker, though more extensive, upper-level outflow- region 
causes a net mass loss in the volume that is reflected in ~° 
the falling surface pressures. As the circulation intensi- 
fies, the level of non-divergence (LND) rises to a higher 
level, and there is a thicker layer of mass influx. 

The horizontal mass flux maxima of the GLAS analysis 
(Fig. 30a) and the ECMWF analysis (Fig. 30b) are consistent 
in magnitude for both surface convergence and upper-level 
@avyvergence. The major difference is that the time of maximum 
Iitateral mass transport in the GLAS analysis lags by 12 h 
the ECMWF analysis. 

The SAT and NOSAT horizontal mass flux fields are shown 


in Fig. 31. As in the analyses, the two-level inflow/outflow 
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regimes are present in both prognoses, and are consistems 
with their respective vertical motion fields. However, as 
the vertical motion fields suggest, the horizontal mass 
fluxes in the prognoses are clearly different from the 
analyses. The NOSAT prognosis (Fig. 31lb) indicates large 
mass outflow (divergence) aloft in the early periods, when 
is consistent with more rapid deepening. During the meee 
explosive phases of oe development, upper-level divergence 
in the analyses exceeds that of the NOSAT, although by Toga 
the intensities of the analysis and NOSAT prognosis are 
Similar. The SAT prognosis (fig. Sla) horizenra wanes. 
transport is similar to that of the NOSAT prognosis, but 
magnitudes of transport are slightly smaller. 

The LND of each prognosis is relatively constant through- 
out development, and is not gradually rising ry: in the 
analysis. This means that in the first 36 h the convergent 
layer was thicker in the prognoses (note position of LND 
at 1318 in each storm). Examination of the LND 1s aygeee 
indicator of each storm's. stage of development. The analysis 
LND increases steadily over the /72-h period and appears 6 
reach maturity at the final time period. The SAT storia 
rises to a maximum after 36 h, levels off and then decreases 
in the final two periods. The NOSAT LND reaches peak |v. tae 
even earlier than the SAT (24 h), although the rise and fag 
are less dramatic than in the other two representations. 


Coincident with the high LND of the NOSAT storm is the maximum 


34 


mass outflow in the upper troposphere. This is most likely 
the result of the NOSAT storm's different relationship to 
the jet streak discussed in Chapter II, and appears to be a 
major factor in explaining early development in the NOSAT 


prognoses. 


fee CONCLUSIONS 

The mass budget study provides explanation for some of 
the differences in development observed in Chapter II between 
the analyses and prognoses. The 12-h time lag in the GLAS 
analysis vertical motion and horizontal mass transport helps 
explain the different analyzed surface pressure in the two 
analyses at 12 GMT 14 January. Aside from this time differ- 
ence, the mass budget results from the two analyses are 
consistent. 

Prognosticmass budgets substantiate overdevelopment 
during the early periods with less intense mass transport 
and vertical motion during the most explosive stages of 
development. The NOSAT prognosis clearly has larger diver- 
gence aloft in the early periods than the SAT prognosis, 
which is consistent with more rapid development. The 
differences in upper-level flow between the prognoses are 
most likely due to the differences in jet stream intensity 
as discussed in Chapter II. 

Static stability throughout all time periods is lower in 
the prognoses than in the analyses, and probably contributed 


to enhanced convective activity and the overdevelopment of 
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both prognoses. This iS particularly true fon the SAg aera 
nosis. Upper-level divergence in the early periods was not 
as strong in the SAT prognosis as in the NOSAT prognosis, 
yet both the SAT and NOSAT prognoses overdeveloped the 
storm. The more rapid destabilization in the early periods 
of the SAT prognosis (Fig. 26) appears to contribute sub- 


Stantially to development. 


2.6 


~ 


Vee Or eer Ty BUDGET ANALYSIS 


A. GENERAL 

The absolute vorticity budget results for the GLAS 
analysis, SAT prognosis and NOSAT prognosis are now pre- 
sented. The Sign convention is positive for processes 
producing vorticity increases in the budget volume (sources), 
and negative for vorticity decreases (sinks). As in the 
previous chapter, all budget figures are for the radius 
10 budget volume. 

The main emphasis of this chapter is a discussion of the 
absolute vorticity time tendency as a measure of storm 
development. The role of various forcing terms such as 
Merrzontal and wacteieal absolute vorticity transport, advec- 
tion, generation and Melee iicn in storm dynamics is ex- 
plored. The horizontal transport terms are separated into 
lateral divergence and lateral advection components to 
investigate their individual contribution to storm develop- 
ment. Vertical transport terms reflect vertical divergence 
and advection from omega values derived from the Kinematic 
method. Table III of Appendix B provides an outline of the 
vorticity budget terms. 

ae in the previous chapter, budget results from four data 
sets are compared. The GLAS analysis is compared to the 


ECMWF analysis as discussed by Calland (1983), and both 
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analyses are used as a basis for comparison of the SAT and 


NOSAT prognoses. 


Bs VORTICITY TENDENG@IESs 

The build-up of absolute vorticity for the cyclone iia 
GLAS analysis over the three-day period is shown in Fig. 
32a. Absolute vorticity increases most notably occur 
between 1300 and 1500. This increase 25 the eesul = so meee 
creased circulation of storm development plus the increased 
planetary vorticity as the storm tracks northward. 

Comparing the GLAS analysis time section with that of 
the ECMWF analysis (Fig. 32b) reveals basically similar fea- 
tures. Slight differences in magnitude and/or time of 
events are observed due to the different resolution of the 
two analyses. 

Vertical time* sections of absolute vorticity for the 
SAT and NOSAT storms appear as Fig. 33. Both prognoses 
indicate an earlier increase in vorticity. In the final 
periods, the SAT prognoses contain relatively higher vorticity 
than the NOSAT prognoses. The decreased vorticity in the 
NOSAT prognoses during the final two periods is consistent 
with the earlier filling of the low in the NOSAT depiction 
of the storm. 

Absolute vorticity time tendency, which is the first term 
in the vorticity budget equation (see Table III), is depicted 
for each of the’ four storms in Figs. 34 and 35. These dis- 


plays are time derivatives of the fields depicted in Figs. 
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32 and 33. The GLAS analysis Vertietty wendency (Fig. 34a) 
is largest in the upper levels for the initial periods, with 
smaller tendencies in the lower layers. The ECMWF analysis 
Mertaicity tendency (Fig. 34b) indicates a minimum in the mid- 
troposphere at 1318, which does not appear in the GLAS analy- 
sis. Calland (1983) attributes the presence of this minimum 
in the ECMWF vorticity tendency to analysis problems, so the 
GLAS analysis tendency is considered to be more accurate 
during this time period. At other time periods the two 
vorticity tendencies are generally consistent. 

Miae SAT vorticity tendency (Fig. 35a) is generally simi- 
lar to the GLAS analysis, except that the key features indi- 
cating certain stages of development clearly occur much 
earlier. SAT vorticity tendency is large in the upper layers 
at 1306, but soon after decreases rapidly ss levels more 
consistent with the analyses. The initially large tendency 
in the jet region is most likely the model's response to the 
weaker winds east of the jet streak. As discussed in Chapter 
II, the magnitude of the SAT jet stream winds at 00 GMT 13 
January was less than the NOSAT prognosis, but by 12 GMT 
13 January the difference had diminished considerably (Figs. 
8a and 9a). Therefore, the model's apparent response to the 
weaker jet winds was to increase vorticity and to strengthen 
the jet. The NOSAT vorticity tendency (Fig. 35b) does not 
indicate the strong initial tendency at the jet level due to 
its different jet winds, but does indicate stronger tendency 


aloft in the middle periods than the SAT prognosis. NOSAT 
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tendency in the lower layers at 1318 is consistent with 


the overdevelopment of the surface low. 


C. LATERAL TRANSPORT OF ABSOLUTE VORTICITY 

The lateral transport of absolute vorticity is one of 
the three major forcing terms in the budget equation (see 
Table III). The lateral transport term can be partitioned 
into components of lateral divergence and advection of 
VOIsEiCi eye 

The GLAS analysis total lateral vorticity transport time 
section (Fig. 36a) indicates strong transport into the 
volume at the upper levels during the initial period, as in 
the ECMWF analysis (Fig. 36b). As circulation becomes well 
established, a transport maxima develops in the surface 
layer that is consistent with horizontal mass transport and 
maximum vertical motion. Coincident with the surface maxima 
is a negative transport maxima aloft. Comparison with the 
ECMWF analysis reveals generally similar features except 
that the GLAS analysis contains a slightly larger inward 
transport maxima in the sunenee layer. Additionally, the 
GLAS analysis exhibits outward transport aloft in the later 
periods, while the ECMWF analysis contains small inward 
iGals e@imumcllO rE: 

Both the SAT (Fig. 37a) and the NOSAT lateral vorticity 
transport (Fig. 37b) more closely resemble the ECMWF analy- 
Sis than the GLAS analysis in terms of structure and magni- 


tude. The major differences between SAT and NOSAT lateral 
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transport are the stronger transport in the upper layers at 
[moo Lor the SAT storm which is consistent with the 
vorticity tendency term. Also, larger transport aloft 


exists in the SAT prognosis in the final periods. 


Pee GATERAL DIVERGENCE COMPONENT 

The Craditional way tO partition lateral transport of 
absolute vorticity is in terms of lateral divergence and 
advection components through the use of vector identities 
(see Table III). The GLAS divergence component (Fig. 38a) 
is composed of a two-layer regime with positive vorticity 
production (convergence) in the surface layer and negative 
vorticity production (divergence) aloft. The surface layer 
convergence maximum occurs at 1418. The ECMWF analysis 
lateral divergence component (Fig. 38b) shows the same two- 
level regime as the GLAS analysis, except that maximum 
surface convergence is less intense and occurs over 12 h 
earlier. Upper-level divergence maxima in the two analyses 
are of similar magnitude and are nearly time consistent. 

Both prognoses lateral divergence (Fig. 39) components 
indicate a rapid increase in divergence aloft in the early 
periods and achieve a maximum that exceeds the analyses. 
The magnitude of the surface convergence maxima in both 
prognoses is slightly less than in the analyses, but the 
maxima are elongated over time. This observation is consis- 
tent with the developmental trends seen in Fig. 1. Both 


prognoses deepen the storm rapidly from the outset and 
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continue at a steady rate until occlusion Seartees ena 
beginning slowly and then undergoing explosive development 

in the middle periods. The NOSAT lateral divergence (Fig. 
39b) indicates larger divergence aloft than in the SAT storm 
(Fig. 39a) from the outset. This stronger surface convergence 
is consistent with the more intense development in the NOSAT 


prognosis. 


E. LATERAL ADVECTION COMPONENT 

The ECMWF vorticity advection component (Fig. 40b) 
exhibits a two-layer regime with negative vorticity advec- 
tion (NVA) below 775 mb and positive vorticity advection 
(PVA) aloft. Calland (1983) attributes the strong NVA 
the surface layer to an asymmetric east-west oriented vorticity 
Maxima coupled to the strong cyclonic flow, which allows the 
advection of smaller vorticity values into the volume. GLAS 
analysis lateral vorticity advection also indicates a two- 
layer regime with NVA in the lower layers and PVA aloft 
(Fig. 40a). The GLAS analysis differs from the ECMWF analy- 
Sis in having relatively weaker PVA aloft and a 24 h lag in 
maximum surface layer NVA. The weaker PVA maxima of the 
GLAS analysis iS consistent with the weaker jet stream winds 
over the surface low at 00 GMT 14 January (Figs. 4a and 10a, 
Calland). 

Strong upper-level PVA during the initial period accounts 
for the large lateral transport of vorticity in the SAT 


prognosis (Fig. 4la). The NOSAT lateral vorticity advection 
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indicates PVA maxima in excess of the SAT prognosis and is 
consistent with the ECMWF analysis. The combination of 
stronger upper-level vorticity divergence and stronger PVA 
is consistent with the more rapid deepening rate of the 


NOSAT prognosis compared to the SAT prognosis. 


F. VERTICAL REDISTRIBUTION 

The third major term in the vorticity budget equation 
(Table III) involves the vertical redistribution of vorticity 
from the lower layers to the middle and upper troposphere. 
The transfer of vorticity in the mid-troposphere intensifies 
mid-level circulation and amplifies the mid-level trough as 
the storm develops. The effect of the vertical transport in 
the vorticity budget is given by the divergence of the verti- 
cal transport. This term can be partitioned into a vertical 
divergence and an advection component. The vertical diver- 
gence component (not shown) is the negative of the lateral 
divergence term discussed in Section D and indicates large 
losses of vorticity in the surface layer and gains aloft. 
The GLAS analysis vertical advection component (Fig. 42a) 
plays a small role in the vertical transport of vorticity. 
A maxima is observed in the mid-troposphere during the later 
periods when both surface layer vorticity and vertical motion 
are large. The ECMWF analysis vertical advection component 
(Fig. 42b) indicates a mid-level maxima similar to the GLAS 


analysis, although it occurs earlier. This is consistent 
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with the differences between the two analyses in the time 
of maximum vertical motion noted earlier. 

The vertical advection components of both prognoses 
(Fig. 43) indicate mid-tropospheric maxima during the later 
periods. The maxima are of similar magnitude to the analyses 
and oce@unm at about the same time as the GLAS analysis (Fig. 
42a). The NOSAT maximm (Fig. 43b) is slightly larger 
than the SAT maximum (Fig. 43a) due to slightly larger verti- 
cal motion in the NOSAT prognosis (Fig. 29). The minor 
contribution of vertical advection to the total budget 
results is consistent with previous studies of DiMego and 


Bosart (1982) and Chen and Bosart (1979). 


C.. SOURGES AND So LLNS 

The final term in the vorticity budget equation (Table 
Milyy) ee the SOLS re input of absolute vyornticiay 
sources and sinks. Sources of vorticity are represented by 
the divergence and tilting terms, and frictional dissipation 
is a vorticity sink. This section will discuss the impact 
of these three components on ecaemeect Ehe four sSceenmn 
representations. 

The generation of vorticity by horizontal divergence is 
mathematically equivalent to the lateral divergence component 
of the lateral transport term discussed in Section C and 
shown in Figs. 38 and 39. Discussion of divergence as a 
source term will not be repeated here. 

The tilting term is the result of vertical vorticity 


generated by the tilting of horizontally oriented vorticity 
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elements by a non-uniform vertical motion field. For this 
term to be sizable, strong vertical shear and a horizontally 
mreving Omega field are required. The tilting terms of the 
GLAS and ECMWF analysis (Fig. 44) are of opposite sign to 
their respective vertical advection terms, which indicates 

a small negative contribution in the mid-troposphere during 
the rapid development phase. Similar features are observed 
maecne SAT and NOSAT prognoses (Fig. 45). 

Frictional dissipation is assumed to occur only in the 
surface boundary layer and is parameterized using a stability 
dependent scheme [Johnson and Downey, 1976]. The frictional 
dissipation of the ECMWF analysis (Fig. 46b) at 1506 GMT 
is of similar magnitude to the low-level maximum in the 
lateral transport of vorticity (Fig. 36b), which suggests 
that frictional effects assume a large role in the dissipa- 
tion of vorticity in the surface layer. The maximum fric- 
tional dissipation in the GLAS analysis (Fig. 46a) is 
slightly higher than that of the ECMWF analysis. Frictional 
dissipation in the SAT and NOSAT prognoses (Fig. 47) are 
also of similar magnitude to their respective lateral trans- 
port terms in the surface layers. Magnitudes of frictional 
dissipation inthe three GLAS fields are in direct proportion 
to the magnitude of lateral mass transport in the surface 


layer and deepening rate in each depiction of the storm. 


fee |6=RESTIDUALS 
The residual in the vorticity budget contains the effects 


of the omitted processes plus the accumulated errors in the 
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calculation of the resolved terms. Enaceuractes in gen 
vertical motion and horizontal wind fields contribute to 

the physical residual components, whereas spatial and temporal 
finite differencing and grid=-point interpolation inaccuracies 
are causes for computational errors [Calland, 1983]. Posi- 
tive residuals indicate that the observed vorticity increases 
at a point exceed those of the computed terms. 

The vorticity budget residual of the GLAS analysis (Fig. 
48a) indicates a sizable vorticity excess in the mid- 
troposphere and a vorticity deficit in the lower tropospijerme 
during the explosive phase of development. Similar resi= 
duals of a smaller magnitude are observed in the ECMWF analy- 
Sis (Fig. 48b). Based on a term-by-term comparison Oia 
two analyses, the larger positive residual at 450 mb 1n@eme 
GLAS analysis is apparently due to an excessive negative 
lateral divergence component. Similarly, the vorticity 
deficit in the lower layers of the GLAS analysis is due to 
an excessive positive contribution from the lateral divergence 
COMpGnenE. 

Residual values for the SAT and NOSAT prognoses (Fig. 49) 
are generally low at all times. The largest residual j6¢ome. 
in the surface layer of the NOSAT prognosis during the early 
periods where observed vorticity tendency exceeds the summem 
the budget terms. This residual is most likely due to the 
different representations of frictional dissipation in the 


model and budget computations. Similarly, different 
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representations of vertical motion in the model and budget 
computations contribute to the residual. As in the analyses, 
the largest residuals in the prognoses occurred during the 
most rapid development phase. For the prognoses, this occurred 
in the early time periods. Development in the prognoses was 
Semcistently strong in the first 48 h, as opposed to modest 
development followed by explosive development in the analyses. 
GLAS analysis residuals resulted from difficulties in the 
resolution of the rotational and divergent components during 
the explosive phase of development. The GLAS prognoses 
residuals resulted primarily from differences between the 
model and budget computational methods for friction and 


vertical velocity. 


m2) 6 CONCLUSIONS 

| Comparison of the vorticity budget Pesarts of the GLAS 
and ECMWF analyses reveals several differences. As in the 
mass budget study, the GLAS analysis surface layer maxima 
of the various vorticity terms consistently lag by 12 h those 
of the ECMWF analysis. Secondly, the GLAS analysis lateral 
divergence of vorticity term is too large in the surface 
layer, and contributes to a negative budget residual. The 
large positive residual in the mid-troposphere of the GLAS 
analysis is most likely due to the 12 h time lag in the 
divergence of the vertical transport. 

Vorticity budgets for the prognoses support earlier and 


more rapid development than in the analyses. The NOSAT 
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prognosis develops the storm more rapidly ebaneence ss 
prognosis primarily due to stronger divergent flow aloft. 
This effect is the direct result of a stronger jet stream 


im the sisiiitialsNosA® tireldss 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A diagnostic verification of the GLAS analysis scheme 
and global forecast model was conducted for a North Pacific 
miso OL ECXtratropical explosive cyclogenesis during 13-16 
January 1979. Storm related mass and vorticity budgets were 
computed for the cyclone using the GLAS FGGE analyses and the 
SAT and NOSAT FGGE prognoses. Results were verified using 
budget results of the ECMWF FGGE analyses as described by 
@arland (1983). 

The coarser grid and more simplified analysis scheme of 
the GLAS analyses compared favorably to the ECMWF analyses. 
Insufficient resolution during the period of maximum storm 
deepening caused the GLAS analysis to under-analyze storm 
intensity at 12 GMT 14 January. At all other times storm 
intensities and positions of the two analyses were very 
close. 

GLAS prognoses initialized with FGGE data forecast the 
storm much more accurately than operational forecasts. The 
Navy operational forecast underdeveloped the storm, whereas 
the GLAS forecasts with FGGE data overdeveloped the storm 
and positioned it more accurately. In this case, the impact 
of the comprehensive FGGE data set on the forecasts was 
positive. 

This case of explosive cyclogenesis was characterized 


by largely zonal mid-tropospheric flow and strong jet streak 
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interaction aloft that induced explosive development. The 
inclusion of underestimated cloud-drift winds in the FGGE 
data weakened the intensity of the jet in the SAT prognesme 
and adversely affected the forecasts. Due to more intense 
jet stream winds, the NOSAT prognosis indicated greater 
surface convergence, greater divergence aloft, greater 
vertical motion and greater PVA than in the SAT prognosis. 
The NOSAT prognosis developed and positioned the storm more 
accurately than the SAT prognosis. The case evaluated in 
this study demonstrates that a larger data base like FGGE 
can improve a forecast. However, the SAT/NOSAT comparison 
suggests that more data may not always be better, and that 
data quality is equally as important as data quantity. 

The approach taken in this thesis was to concentrate 
on a single system and perform an in-depth diagnostic study 
and describe the processes involved in explosive cyclogene- 
Sis. The conclusions reached in this thesis are for this 
Single storm. It is recommended that future research examine 
other extratropical maritime explosive cyclones and specifically 
cases for which both SAT and NOSAT prognoses are available. 

Without the FGGE, a study of this type would not have 
been possible. Data gaps have always made forecasting over 
ocean areas difficult. Satellite data in recent years have 
helped greatly in filling data gaps, and imagery has unques- 
tionably aided the forecaster. It is not known if the errant 


jet stream cloud-drift winds seen in this study are unique 


a0 


to this particular storm. It 1s recommended that the 
procedure for evaluating cloud-drift winds prior to data 


assimilation be reviewed. 


Se 


APPENDIX A 


FGGE DATA 


The data collected during the first special observation 
period (SOP-1) of the Global Weather Experiment represented 
the most intensive effort ever attempted in atmospheric and 
oceanographic data collection. The experiment clearly 
demonstrated that state-of-the-art technology coupled with 
outstanding resource management can produce the higher quality 
data sets needed for the improved performance of numerical 
prediction models. 

The GLAS objective analysis procedure is characterized 
by Baker (1983) as a modified Cressman scheme based on a 
method developed by Bergthorsson and Doos (1955). The GLAS 
assimilation/forecast model consists of nine vertical levels 
in sigma coordinates with a non-staggered horizontal resolu- 
tion of 4 degrees latitude by 5 degrees longitude. All 
horizontal differences are computed with fourth order 
accuracy. 

Substantial preparation of the FGGE level II-b data was 
performed by GLAS prior to implementing the objective analy- 
Sis scheme. Of the changes made to the level II-b data, 
the most significant to the study of the North Pacific ™Seonm 
was the correction of the Japanese cloud-drift winds. The 


Japanese upper tropospheric cloud-drift winds were considerably 
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Mice in Magnitude tham surrounding rawinsonde and air- 
craft reports. This error was attributed to erroneous 
altitude assignment [McPherson, 1984]. Altitude errors were 
likewise observed in the Japanese lower tropospheric cloud- 
drift winds. The method used was to assign the wind alti- 
tude to the low cumulus cloud-top altitude as determined by 
brightness temperature. Studies determined that low- 

level cloud-Grift winds are m@st accuarate at the cloud base, 
and not the cloud top [McPherson, 1984]. Consequently, all low 
level cloud-drift winds were reassigned by GLAS to the 

900 mb level (the average cloud-base level over the ocean). 
Altitudes of cloud-drift winds between 400 mb and the tropo- 
Pause were reassigned based on cloud-top brightness tempera- 
ture [Baker, 1983]. 

The GLAS analysis first-guess is provided by the 6-h 
model prognosis at the model sigma levels and must be inter- 
eeeaecd tO the mandatory pressure levels. Observed data are 
then ingested and compared to the first-guess values. Data 
differing from the first-guess by a specified amount are 
flagged as suspect and then checked against surrounding data. 
If acceptance criteria are still not met, the data are 
rejected. The "good" data are then interpolated to the 
4 degree by 5 degree grid using the Cressman scheme, which 
has been modified to consider data quality and variable data 
density. The u and v wind components, heights, and relative 


humidity are analyzed at the mandatory pressure levels, and 


5.5 


the surface pressure and surface temperature are reduced to 


sea level and analyzed. 
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APPENDIX B 


TABLES 


TABLE [I 


Generalized QLD Budget Equation in Isobaric Coordinates 
(After Wash, 1978) 


Re os em a > 
F = i (= eoencto cd de (—dp) 
1000 mb 0 0 J 


where F is the volume integral of the desired budget property 
fee 


The budget equation is 
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where the lateral transport is 
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and the vertical redistribution is 
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The source/sink term is 
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QLD Mass Budget Equations in Isobaric Coordinates 
(After Wash, 1978) 


The Definition 
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TABLE rt 


QLD Circulation Budget Equations in Isobarie Coordinage=. 
(After Wash, 1978) 
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where Co is the absolute circulation and 5 is absolute 


vorticity. The budget equation is 


—_ = LT(c_) + DVT (z_) + S(t.) 


where the lateral transport is 
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and the divergence of the vertical transport is 
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The source/sink term is 
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TAB eat 7 (CONT. ) 


Section B 


The partitioned form of the vorticity budget equation is 





S(C_) 
= - 
FE er AC) + DVT (c_) S (c_) 
vertical vertical 
divergence advection 
Memizontal horizontal divergence tilting frictional 
divergence advection term term dissipation 


The above partitions make use of Stokes' theorem 
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Synoptic Fields for GLAS Analysis at 0000 GMT 

12 January 1979. (A) 250 mb winds/Isotachs 
(m/s). (B) 500 mb Absolute Verticity (DASHED) 
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(C) Sea Level Pressure (SOLID) in mb and 1000-500 
mb Thickness (DASHED) in gpm. . 
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Figure 26. Stability Index. Values in Degrees Kelvin/ 
100 mb. Time 1300 Refers to 00 GMT 13 January. 
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Figure 28. Analyses Kinematic Vertical Velocities 


(A) GLAS (B) ECMWF. Units are in mb/1000-sec. 
Time 1300 Refers to 00 GMT 13 January. 
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Figure 32. Analyses Area-Averaged Absolute Vorticity | 
: Vertical Time Sections (A) GLAS (B) ECMWF. ‘Time 
1300 Refers to 00 GMT 13 January. Units are 
ie A= 57SEC. | 
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Figure 35. Prognoses Absq@lute Veoreicity Vertical Tie 
Tendency Sections (A) SAT (B) NOSAT. As in 
Figure 34. 
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Figure 36. Analyses Lateral Vorticity Transport 
(apmeGeAGe(s) ECMWE.. AS in Figure 34. 
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Figure 39. Prognosés Laperameyortirelty pn vergenc= 
(A) SAT (B) NOSAT. As in Figure 34. 
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Ficure 40. Analyses Lateral Vorticity Advection 


Pemereoe(8) ECMWP. 9 As in Figure 34. 
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Figure 41. Prognoses Lateral Vorticity Adcvection 
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Figure 42. 
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Figure 43. 
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Figure 44. Analyses Tiltine Term (A) GLAS (8) ECMWF. 
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49.  Prognoses Budget Residuals (A) SAT (B) NOSAT. 
As in Figure 34. Positive Values Indicate 
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